Effect of environmental coupling on tunneling of quasiparticles in Josephson junctions 
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We study quasiparticle tunneling in Josephson tunnel junctions embedded in an electromagnetic 
environment. We identify tunneling processes that transfer electrical charge and couple to the 
environment in a way similar to that of normal electrons, and processes that mix electrons and 
holes and are thus creating charge superpositions. The latter are sensitive to the phase difference 
between the superconductor and are thus limited by phase diffusion even at zero temperature. We 
show that that term is suppressed in many environments, thus leading to lower quasiparticle decay 
rates and thus better superconductor qubit coherence than previously expected. 
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The physics of micro and nanoscale Josephson junc- 
tions is a paradigmatic application of macroscopic dissi- 
pative quantum mechanics of open systemiP^. This, on 
the one hand, makes them ideal test-beds for that theory, 
with largely tunable parameters^]. On the other hand, 
Josephson junctions have an ample range of applications 
in sensing^, amplification-^, metrology® and quantum in- 
formation processings^. p or these applications, it is im- 
perative to thoroughly understand the dissipative quan- 
tum physics of Josephson Junctions in order to achieve 
optimal device performance. 

Most older work on this topic has focused on the trans- 
port properties of Josephson junctions, their current- 
voltage and noise characteristics. The 'P(E) '-theory of 
treating environmental fluctuations^ 3 -' has emerged as a 
powerful tool in the case of small junctions when (Joseph- 
son or quasiparticle) tunneling can be treated perturba- 
tively. 

For applications in quantum computing, another ques- 
tion in this framework has emerged 14 * 15 ! Rather than 
computing current- voltage characteristic, the total quasi- 
particle transition rate is the quantity of interest. This is 
important because quasiparticle transitions in any direc- 
tion are highly detrimental to qubit coherence 16 . Refer- 
ences |l#fT5] study this problem in an approach that is 
perturbative in the coupling to the environment as mea- 
sured by the environmental impedance in units of the 
quantum resistance Rq — h 2 /e ~ 25.5 kfl. In this paper 
we are going to generalize that study to arbitrary system- 
environment interaction strength. We are going to show 
that some of the processes mix the electron-like and hole- 
like branches of the quasiparticle spectrum. These are 
generally phase-dependent but easily suppressed by the 
environment. These observations lead to slower relax- 
ation than expected in previous work. 

We formulate the mathematical description of quasi- 
particle tunneling in an environment. We start from the 
Hamiltonian 
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FIG. 1: Quasiparticle processes labeled by their dispersion 
connected to an external circuit. By linearization, the Joseph- 
son circuit is replaced by an effective impedance. 



Here, we have the BCS Hamiltonian d 19 ' 2 Q in mean field 
form describing the two electrodes 
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where = is the kinetic energy of the electron, /Zj is 
the chemical potential in superconductor i = 1,2, being 
Hi = Ep and \i2 = Ep + eV where V is the applied 
voltage. We keep Aj = |A| e 1 ^ 1 / 2 complex in order to 
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allow for a phase bias. The tunneling Hamiltonian is 

= ( T kA,a.l d l^ + T fci4,<x,2^,l) ( 3 ) 
kl 

Before discussing the electromagnetic environment, we 
diagonalize the BCS Hamiltonians, eq. ^ through the 
Bogoliubov transformation: 

Cfct = «fc7fct + v kll kV 
C-ki = ~ v *k%t + U kl- H - 
Here, the BCS coherence factors are defined 
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where we introduced the quasiparticle energy E^i = 
\J (£k,i ~ Mi) 2 + l^l 2 - This allows us to rewrite the tun- 
neling Hamiltonian, cq. [3] in the following form 



x (tLi7^2 + 7l ka2 l-iai) + H T2 (6) 



where we introduced the phase difference (j) — <j>\ — <j>2- 
The term Ht 2 contains operators that change the num- 
ber of quasiparticles by two, i.e., contain terms of the 
structure 77 and 7T7T hence changing the total number 
of quasiparticles in the setup, which do not contribute to 
the quasiparticle rate -these terms contribute to Joseph- 
son and Andreev processes. We now want to evaluate the 
Fermi's golden rule rate for a transition that transfers a 
quasiparticle from the electrode 1 to the electrode 2. The 
relevant matrix element is 
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We now use that for a nonmagnetic barrier, T) 
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T*u_, can be chosen real. In the Fermi's golden rule 
transition rate, we need the absolute square 
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where u(E,E') = \ui(E)u 2 (E')\ and v = M£> 2 (£")l 
with Ui/ 2 (E) and Vi/ 2 (E) given by the BCS coherence 
formulae, eqs. Q and §5§ Thus, the transition prob- 
ability contains an interference term that is sensitive to 



the phase bias of the junction. This term describes a pro- 
cess that, even though it transfers a genuine quasiparticle 
from electrode 1 to electrode 2, it actually consists of a 
superposition of electron and hole transfer, i.e., it is not 
diagonal in charge space, see the diagrammatic represen- 
tation in fig. [2] We will see later on how this is important 
in the sensitivity to environment. This term bears anal- 
ogy to the famous cos <fi term in the Josephson effeciPStHl 
Refs. |14I15| have studied the unbiased, = limit of 
this formula. 
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FIG. 2: Diagrammatic representation of the different contri- 
butions to the quasiparticle tunneling matrix elements. Left: 
Processes closed in quasiparticle and in charge space, right: 
quasiparticle process that creates charge superpositions. The 
typical state |n, m) represents n particles in the lead 1 and m 
particles in the lead 2. 

In order to capture the influence of the environment, 
we apply the ideas of P(i?)-Theory^l There, the envi- 
ronmental Hamiltonian is described by an oscillator bath 
-ffbath = X)n u, n&n<hi and couple the oscillators linearly 
to our quasiparticle system H env - C — N\e5V , where we 
assume that only the first reservoir fluctuates - this cor- 
responds to a specific choice of gauge. The total number 
and voltage operators are 



(7) 



We can now work out the total tunneling rate by sum- 
ming over momentum states at constant energy from the 
initial state \i) = \R)\E) to the final state |/) = \R')\E') 
and considering that the matrix element of environmen- 
tal Hamiltonian (E\Tk q r yq'l'y 2 . cr \E') thermal distribution 
of quasiparticles. The transition rate formally can be 
written 
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where Di is the reduced density of states and fa is the cor- 
responding distribution functions on the lead i. We con- 
sider here that the nonequilibrium quasiparticle tunnel- 
ing keeps the system in thermal and chemical-potential 
equilibrium. The influence of the external circuitry on 
quasiparticle tunneling is encoded in the dressed vertex 
P tot (e), the probability density for exchanging energy e 
with the environment. From the Fermi's golden rule the 
tunneling transition rate can be written in more detail as 
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fi(V) = -2— / dEdE'D 1 (E)D 2 (E')f 1 {E)(l- f 2 {E'))J2 \(R'\H T \R)\ 2 P ? (R)5(E + eV - E' + E R - E R ,) 



Pp(R) = (R\pp\R) for Pfs = e -P H °™ /X)e-^»" and the 
arrow on T indicates the direction of tunneling from the 
lead 1 to 2. Rn is the tunnel resistance of the junction in 
the normal state. By expanding the Dirac delta as a tem- 
poral integral 6(x) — (2irh) 1 dte~^ xt and writing 

the operator (R'\e~ E ^' t H T e ERt \R) = (R'\H T (t)\R) the 
time evolution of the Hamiltonian makes the rate to be- 
come proportional to ((Hj. (t) H!^ (0))), where ((• • • )) = 
Sfl R' (R' \ " ' Pp\R)- Substituting the tunneling Hamil- 
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tonian into this formula one gets ((Ht (t) Ht (0))) = 
(((ue^W/ 2 - ue -^W/ 2 ) ( ue -^( )/ 2 - ve**^ 2 ))). 
Taylor expanding the exponents and reordering it 
is easy to check (e±**(*) e =F«*(o)) _ e <(*(0-*(0))*(0)) 
and ( e ± #(*) e ±i ^( )) ~ e -(Wt)+^0W(o)>, Sus- 
bstituting this one can get (Hi, (t) H^ (0)) = 

U u 2 + y 2\ e (0(t)0(O)>/4 _ 2uue -(0(t)«i(O))/4\ e -(^(O)0(O)>/4_ 



This makes the transition rate to become 
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where u(E,E') = \ Ul (E)u 2 (E')\ and « = \vi(E)v 2 (E')\ 
with u 1 / 2 (E) and Vi/ 2 (E) given by the BCS formulae, 
eqs. @ and ([5|. 

The vertex probability is defined as the probability 
density for exchanging energy e between system and the 
environment and is defined formally as 



On the other hand, J a (t) = ((6(j>(t) + 5<f){0)) 6 (f)(0)) /4 = 
J s (t) + (<50 2 (O)) /2 captures processes that mix electrons 
and holes. 

Other than J s , this correlation function does not van- 
ish at t — 0, rather, J a (0) — (<50 2 (O)) /2 captures zero- 
point fluctuation. Consequently, 
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one can correctly ex- 
■ 2 ) e- J ^ -2uve~ J ^ t \ 



By comparing eq. Q and Q < 
pect the definition e J &> = (u 2 + v 2 
where e Js(t) = ( e ^( t )/ 2 e~^ (0)/2 ) ~ e -<(0(t)-0(o)Wo))/4 
and e J ° (t) = ( e *( t )/ 2 e l0(o)/2 ) - e ^ ( * (t)+ ^ o))0(o)>/4 . 
From these equalities one can simplify the definitions 
to be J s (t) = ((S(f)(t) - S(f>(0))5(f>(0)) /4 and J a (t) = 
— ((S(f)(t) + 5(f)(0))S(f)(0)) /4, where the phase factors are 
given by 



(11) 



S(f>(t) = - / dt' 5V(t'). 



Note that in the definitions we considered that the mean 
value of the phase is not physically relevant as long as the 
effect of environment is summarized into the presence of 
phase fluctuations about this mean value. The symmetric 
combination J s (t) captures all-electron and all-hole pro- 
cesses and the corresponding vertex coincides with the 
all-electron P(E) known from Ref|3] 



P S (E) = P(E) 



dte~ iEt e 



Js(t) 



(12) 



Pg(E) 
Ps(E) 



(^ 2 (0)\ _ 



-Et 



(i<f> 2 (0)) 



tfdte J W+i 



EL 



(13) 

P S (E), P a (E) are not normalized as they are density 
probabilities. We can understand this as follows: Even 
though electron-hole mixing processes are diagonal in 
the Fock space of quasiparticles, they are off-diagonal 
in charge space. The electromagnetic noise couples to 
these electrical charges. The environmental modes dress 
the charge in an attempt to measure charge and local- 
ize it. This has an analogy: Creation of a superposition 
of charge states occurs, following the charge-flux uncer- 
tainty relation, whenever the phase across a Josephson 
junction is becoming localized, i.e., when a Josephson 
junction is behaving classically. Thus, on the one hand, 
the scaling factor in eq. (13 1, directly measures the de- 
gree of charge localization. On the other hand, research 
aiming at maximum supcrcurrent (hence maximally lo- 
calized phase and maximally extended charge) in small 
Josephson junctions in a dissipative environment finds 
exp(-(</>(0) 2 )/2) to be the relevant reduction factor^Ml. 

In order to compute J s / a , we can show, by introduc- 
ing an oscillator bath model and matching its proper- 
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ties to the fluctuation-dissipation theorem as it describes 
Johnson-Nyquist noise. Applying the standard proce- 
dures used in P(i?)-theory, we find 
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Notably, this integral seems to diverge at its infrared end, 
where it formally appears to be ~ J ^j-Z c h (0) at T = 
and oc T J ^ReZ eS (0) at T > 0. In the combination 
J s (t) = S(t) — 5(0) this divergence is regulated, but the 
same argument does not apply to J a (t) = S(t) + S(0). 
Thus, if Z c s(0) 0, eq. (13) enforces P a — 0, meaning 



that all electron-hole mixed processes would disappear. 
We will see later that this does not occur in standard 
physical environments. 

The final expression for the probability density 
P{E, E') is thus 



P^{E,E>) = J ^Le-^V-z'K 
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In order to evaluate the typical P s / a (E) and the 
resulting tunneling rates, we need to project likely 
models of the electromagnetic environment encoded in 
Z c g. Fol lowing the RCSJ models and its microscopic 
analogue j 1 ! 2 ! 27 ! 28 !, we consider the quasiparticle channel 
to be put in parallel to the supercurrent, Josephson chan- 
nel, the junction capacitance, and an external impedance 
Z(u), see Fig. ([I]). We model the linear impedance 
of the Josephson channel by its Josephson inductance 
Lj = $o/(27r/ c ). The total effective impedance of this 
parallel setup leads to 
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is the junction's plasma frequency. 



In the case of a simple resistor or lossless transmission 



line, Z{iS) — R and we find 
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Here, we can observe that at low frequencies, Re.Z e ff oc 
u! 2 . Physically, the reason for this is that the noise from 
the resistor gets shunted through the Josephson junction 
at low frequencies and does not affect the quasiparticle 
channel. Thus, although P a will be normalized to a value 
smaller than unity, the scenario of it scaling all the way to 
zero does not occur. This is consistent with the fact that 
phase qubits do show a zero-voltage current. The same 
general conclusion holds for other physical environments 
tested. 

Given the high quality of qubit junctions, it is appro- 
priate to start from the limit of an infinite quality factor, 
R — > oo. In that case, the effective impedance reduces to 



Z c s = ^ [S (u) - uj p ) + 5 (w + uj p )] . (18) 



Introducing the dimensionless parameter p — Tr^p, 

where Z = \J ^4- is the wave impedance of 
the junction, and given the definition of Z e tt — 



ir/2C (5(w — uip) + S(lu + cu p )) and that u> p = 1/y/LjC 
and Lj = $>q/2ttI c we can find that 5(0)/2 = p = 
ttZq/ Rk- For the total environmental probability we use 
the result of RefP 



P ( E ) = Pk(p,u P T)6(u - kcjp) 
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which has the form of a series of sidebands, corresponding 
to the emission / absorption of k photons to/from the 
plasma mode. The weights in the low temperature limit 



(16) of ksT <§C Hujp is Poissonian pk (p,oj p ) 



- p p k /k\. This 



coincides with that of^in the limit of p — > 0. The total 
quasiparticle tunneling rate hence reads 
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For large capacitance junctions where the phase thus one can get the approximate result at low phase 
is localized the exponential can be Taylor expanded, limit by the approximation (u 2 + v 2 ) exp [S(t)/2] — 
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2uvexp [S(t)/2] ~ (it - vf + (u + vf S(t)/2. This ap- 
proximate limit result is consistent with the results taken 
from RefPS' except that ours show that there is a prefac- 
tor exp(— (S<j>(0) 2 ) /2) in the rate that is missing in previ- 
ous works that can reduce the quasiparticle environmen- 
tally mediated transition rate. 

Note that for the special case of n = 1 we find 
rbaro,ie _p + rdr,i e ^ where rb are ,n corresponds to the 
bare quasiparticle tunneling rate 
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under a bias voltage controlled by the Josephson-type 
relation V — nhu>/e and 
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RC \nu J D(E F )A 
where D(Ep) is the density of states in the leads at the 



Fermi energy and n qp the density of quasiparticles in the 
sample. 

Let us now compute the rate. For this to be done 
we need to know the explicit energy dependence of the 
transition rate. Using eq. Q and ^ one can show that 
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Substituting the density of states and the coherence 
factors in the transition rate we arrive at the complete un- 
perturbed formula for nonequilibrium quasiparticle tun- 
neling 
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This is one of the main results of this paper based on 
which one can calculate the rate for any subsystem. As 
one can see the asymmetric term P a (E) is reduced, ac- 
cording to eq. [l3] by a factor of e~ 2p . This leads to a 
difference in O(p) that occurs in our approach account- 
ing for zero-point fluctuations in the impedance as well as 
phase fluctuations. For phase qubits, whose impedance 
is engineered to be Z ~ 50S7 this correction does not 
qualitatively change the physics, however, improves the 
quality of fits to the data. In other type of qubits with 
higher impedance these corrections will be crucial. A 
qubit relaxation/excitation experiments involving non- 
equilibrium quasiparticles^, will attempt to relax the 
qubit energy splitting which is ~ uj p into quasiparticles, 
hence, only the main band is relevant. 

Overall normalization of P(E) means that the appear- 
ance of sidebands generally lowers the main resonance. 
In order to study other types of qubits, we can easily 
show that 



E c 
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For a transmon, Ej/E c rs 30 and for a flux qubit ~ 
50, making this correction large enough to be visible in 
a reduction of the quasiparticle rate whereas the cos <j>- 
term will not be suppressed altogether. More explicitly, 
for a transmon we get a deviation of around 7% from 
the perturbative approximation (1 — e~ 2p sa 2p) and for 



a flux qubit we get a deviation of around 5%. For a 
traditional charge qubit, p is large, leading to a different 
regime. FroirpSl Ej/E c ss 0.35 in a typical charge qubit 
case leading to a large deviation of around 40% from the 
perturbative approximation. 

We are now studying the opposite case, an environmen- 
tal mode that is overdamped by an external impedance. 
Overdamping means the width of the resonance in Z e ff, 
eq. 17 is larger than its frequency, i.e. uj p RC <C 1 or 
equivalently R -c Zq. In that case, to lowest order in the 
small parameter r = ^ is, the poles of Z e s, eq. |l7]are at 
uj = ±17 and at cj = ±ijr 2 where 7 = 1/RC. Note that 
unlike other work for overdamped oscillators^We keep 
the next-to leading order which keeps the second set of 
poles to zero, rendering the quasiparticle current nonzero 
in view of eq. (15 1. This will have an important conse- 
quence later on. 
rewrite eq. ( 17 1 as 



We can now, in the overdamped case, 
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This can be read as the difference between two Ohmic 
environments with a Prude c utofPSl. This allows us to 
rely on known results^*- 5 * 13 * 30 ! for computing J(t). Here, 
we focus on the zero-temperature regime. We find for the 
zero-point term S(0) = —-jpt logr. 

Now, interestingly, this wSl be a large positive term for 
r — > which can be achieved by going to small u) p , i.e. 
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the uv -term in eq. ( 15 1 will be suppressed by a factor 
r iR/R K t a very small value. 

It is known 13 is that for = , , R P, 2 we have 



"E 1 (-w T)-eP aT E 1 (w Tj\ (26) 



and that for long times, loqt — > oo this reduces to Js(t) = 
(logw c t +7 e + if) where 7 e ~ 0.5772 is the Euler- 
Mascheroni constant. For short times, lu^t <C 1, we find 
J(t) = iiraujQTrt. We can now find P(E) in three regimes. 
For E <C 7,7r 2 , both terms in eq. 25 should be treated 
in the long-time limit. In lowest order in r we find 



J(t) 



4R 
Rk 



logr 



(27) 



i.e., no time-dependence, in agreement with the fact that 
the env iron ment is super-Ohmic at these low frequencies, 
" This leads to P„(E) = r 4 PS(E). 
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see eq. 

In the intermediate regime, jr 2 <C E <C 7, we can 
combine short and long-time limit as 



J(t) 



Rp 
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log jt + 7 e + i — — iirjr t 
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This leads to an Ohmic P(E) that is energetically 
shifted by SE = jf^v^r 2 = This results in 



P(E) = 



T(2R K /R) E V R E, 
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where we introduced the charging energy E c = ^=?. Fi- 
nally, at large energies hence entirely short times and 
large energies, E 3> 7, we can approximate 



J(t) ~ -i 



(30) 



hence leading to a simple capacitive contribution 
from the junction's charging energy and P{E,E') — 
e -s(o) (( u 2 + v 2) e s(t) _ 2uve- s W) S (E' — E — e 2 /2C). 

In conclusion, we have developed the theory of quasi- 
particle tunneling for superconducting tunnel junctions 
in an arbitrary linear dissipative environment. We 
worked out the unperturbed tunneling rate of nonequilib- 
rium quasiparticles in the junction. In the perturbation 
regime P(E) governs processes of an electron or hole tun- 
neling, which are represented by the diagonal transition 
matrix elements in charge space. We showed that pro- 
cesses that create charge superpositions are additionally 
supressed by zero-point fluctuations of the phase. For 
the case of low damping, quasiparticle tunneling is ex- 
changing energy with the plasma mode in the form of a 
sequence of sidebands. 
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